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Weak topological insulators have an even number of Dirac cones in their surface spectrum and 
are thought to be unstable to disorder, which leads to an insulating surface. Here we argue that the 
presence of disorder alone will not localize the surface states, rather; the presence of a time-reversal 
symmetric mass term is required for localization. Through numerical simulations, we show that in 
the absence of the mass term the surface always flow to a stable metallic phase and the conductivity 
obeys a one-parameter scaling relation, just as in the case of a strong topological insulator surface. 
With the inclusion of the mass, the transport properties of the surface of a weak topological insulator 
follow a two-parameter scaling form. 

PACS numbers: 73.20.-r, 71.23.-k, 72.15.Rn 



The conventional band theory categorizes crystals as 
metals, semiconductors, or insulators depending on the 
size of their band gap. Over the last few years it has been 
realized that this categorization overlooks the fundamen- 
tal fact that not all insulators are equivalent. Topological 
insulators are characterized by nontrivial band topology 
leading to gapless metallic surface states which are ro- 
bust to disorder that preserves time-reversal symmetry 
(TRS) [lj. In two dimensions (2D), the quantum spin 
Hall (QSH) insulator possesses a pair of counterpropagat- 
ing edge modes which are protected from backscattering 
by TRS. In three dimensions (3D), topological insulators 
are classified as either strong (STI) or weak topological 
insulators (WTI). The surfaces of STIs have an odd num- 
ber of 2D Dirac fermions and have garnered much of the 
attention, as TRS disorder cannot localize the surface 
states unless it is strong enough to move states across 
the bulk energy gap. In contrast, the WTIs have an even 
number of Dirac fermions and are believed to be unstable 
to disorder 

This belief stems partially from comparisons with 
graphene. Superficially, WTIs and graphene are simi- 
lar in that their low energy electronic properties are de- 
scribed by an even number of Dirac fermions @] . While 
both systems have TRS implemented by an antiunitary 
time-reversal operator O, they differ fundamentally in 
that 6 2 = +1 for graphene from SU(2) spin symme- 
try while for a WTI 6 2 = — 1 due to the presence of 
strong spin-orbit coupling. This places graphene in the 
orthogonal (AI) symmetry class, while WTIs belong to 
the symplectic (All) class in the Altland-Zirnbauer clas- 
sification 0] ■ The consequences of the minus sign are pro- 
found. The first quantum correction to the Drude con- 
ductivity is determined by interference of time-reversal- 
symmetric paths. In the orthogonal class, this interfer- 
ence is constructive (weak localization) and eventually 
leads to localization of all single particle states. In con- 
trast, in the symplectic class, the interference is destruc- 



tive (weak antilocalization) giving rise to an enhance- 
ment of the conductivity and a stable symplectic metal 
phase [1| . Hence, the metallic phase of graphene is un- 
stable to disorder coupling the Dirac fermions @ but is 
stable in WTIs. 

An STI is also in the symplectic class. With an odd 
number of Dirac fermions on its surface, it always flows 
into the symplectic metal 0, HI , reflecting the presence 
of a topological term in the effective field theory (non- 
linear sigma model) describing diffusion 0, [TcJ. This 
topological term is absent in the same description of a 
WTI, suggesting that localization should occur. In con- 
ventional semiconductors with spin-orbit coupling, this 
leads to a metal-insulator transition at a critical conduc- 
tivity a c w 1.42 e 2 /h 
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It is the purpose of this work to explore the precise 
conditions under which a WTI undergoes localization. 
One reason that this is a pressing question is the follow- 
ing argument 12j . If one considers obtaining a WTI by 



stacking 2D layers in the QSH phase, a surface parallel 
to the stacking direction would consist of pairs of one- 
dimensional (ID) countcrpropagating helical modes. The 
number of such modes taking part in transport can be 
even or odd depending on the number of layers. However, 
an odd number of ID modes in the symplectic class nec- 
essarily leads to the presence of a perfectly transmitted 
mode and thus a minimum conductance of e 2 /h [l3, 14 1. 
While this argument is one-dimensional in nature as the 
sample thickness is constant, it suggests that a WTI can 
under certain conditions avoid localization. In the ex- 
tended two-dimensional surface, the meaning of this par- 
ity effect is unclear, raising the question: What is the 
scaling behavior of the conductivity in disordered WTIs? 

In this Letter, we demonstrate, by numerical simula- 
tions, that the scaling flow depends on the presence or 
absence of a specific time-reversal-symmetric mass, to be 
defined below. In the presence of this mass, a gap opens 
up in the spectrum which can lead to localization. Dis- 
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Disorder structure 



Disorder type 



Notation 



Ko ■ t x 
V yx ■ r y a x 
V ■ T v n v 

Vyy T a 

V yz ■ r y a* 
Voo ■ t 



scalar potential (2xAII) 
gauge potential (2xAIII) 
gauge potential (2xAIII) 

mass (2xD) m = {V yz ) 

scalar potential (2xAII) 
scalar potential (2xAII) /j, — —(Voo) 



TABLE I. List of time-reversal invariant disorder terms on 
the surface of a WTI with two Dirac cones. If only one of the 
disorder structures is present in the system, the type indicates 
the disorder class of the system and the effect of the disorder. 
For example, with only V yz (r)T y a z , the system breaks up into 
two systems, each identical to a Dirac cone with random mass 
in class D. Hence, multiple disorder structures are required for 
the system to be class AIL 



order can still drive the system into a metallic phase, 
realizing a metal-insulator transition at a critical value 
of conductivity consistent with what is observed in con- 
ventional semiconductors. In contrast, in the absence 
of this mass the system always flows into the symplcc- 
tic metal. We demonstrate that this flow follows one- 
parameter scaling with a positive beta function, just as 
in the case of an STI 0, [§] . The phase diagram emerg- 
ing from these observations (cf. Fig. [T|) suggests that 
one-parameter scaling is not realized throughout, as one 
might expect from the minimal nonlinear sigma model 
description. Instead, we present data supporting two- 
parameter scaling, the effective field theory of which re- 
mains unknown. 

Hamiltonian and disorder structure. In the following, 
we specialize to the case of a WTI with two Dirac cones, 
for which the low energy electronic properties are de- 
scribed by the Hamiltonian 15 1 



H = hv D T°{a x k x + a v k y ) + V(r), 



(1) 



where t° = ct = 1 is the identity, t x ' v ' z and (j x ' v,z are 
the Pauli matrices in valley and spin space, respectively. 
H is invariant under the time-reversal = ia v JC, where 
K, is the complex conjugation operator. The Dirac veloc- 
ity Vd (taken isotropic for simplicity) and H are set to 1 
henceforth. The disorder potential is written 



V(r)="£v a p(r)T a 



(2) 



at) 



with V a /3 (r) a scalar potential and a, ft E {0, x, y, z}. The 
six terms respecting time-reversal, listed in Table HI are 
independently distributed with correlation 



{5V aP {v) 5V a p(v')) = g aP K{r - r') 



(3) 



where Jd 2 rK(r) = 1. The two-terminal conductiv- 
ity g of a system of size L is obtained numerically 
by adapting the transfer matrix method of Ref. to 




FIG. 1. The phase diagram of the Hamiltonian fT} as a func- 
tion of mass m and disorder strength g a p — g. The solid 
line marks the metal-insulator transition at /i = 0, whereas 
the dashed line marks the transition at finite fi. At the clean 
Dirac point (g = fi — 0), there is a topological phase tran- 
sition between the two types of insulators. With increasing 
disorder or chemical potential fj,, a metallic phase appears 
separating the two topological sectors. 



the current problem. (The width W is taken large 
enough that the conductivity is independent of the ratio 
W/L.) Each disorder term is Gaussian correlated with 
K(r) = cxp(— ?- 2 /2£ 2 )/(27r£ 2 ). We also take the averages 
{V a p(r)) = 0, except for V yz and Voo as explained below. 

It is useful to first analyze the system in the clean 
case, where V a /s are constants. (Vbo) acts as the chemical 
potential /j, which shifts the energy spectrum trivially. 
T y o~ z anticommutcs with all the other potentials (except 
1) as well as the kinetic term a ■ k; the presence of this 
term always gaps the system, and hence we refer to m = 
(V yz ) as the "mass." The energy spectrum of the system 
is given by 



- 2 i y2 1 , y2 ,y2 ,y 



± 2 ^ (V X 2 + V 2 2 0)fc 2 + (V yx k x + Vyykyf + HI*, (4) 



with minima at k 2 = V x0 



V 2 +V 



i V^o and k x /ky 



Vyx/Vyy, in which case we have (E — fi) 2 = rn 2 . There- 
fore, the energy gap is 2\m\ and the system is insulating 
when \m\ > The cases m > and m < — \fi\ corre- 
spond to the two topological sectors in the 2D All class, 
i.e., the trivial and QSH insulator [16j . The intermediate 
metallic region — \fx\ < in < \fj.\ separates the two phases. 

In the presence of disorder, a similar description ap- 
plies - by varying to, one can take the system between 
the two insulating phases. As conjugation by t x flips the 
sign of to, a conducting state should be realized at to = 0. 
Because to the stability of the symplcctic metal, one does 
not expect genericallya direct transition between the in- 
sulating phases p7Hl9j . The resulting phase diagram is 
shown in Fig.Q] The shape of the phase diagram around 
the clean Dirac point g = m = fx = 0is consistent with 
the renormalization group flow of the coupling parame- 
ters g a/ 3, to, and [i away from that point (c/ |App. "A"| . At 
a finite chemical potential, there is a range of mass val- 



ues TO 



< 



\x\ where the system undergoes two transitions 
with increasing disorder strength (dashed line in Fig. [1]) . 
A similar mass term can be defined for an arbitrary even 
number of Dirac cones; thus, the phase diagram in Fig.Q] 
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FIG. 2. Demonstration of one-parameter scaling at m = 0. 
Conductivity as a function of system size for various parame- 
ters all collapsed (by shifting the raw data horizontally) onto 
one scaling curve. At large a, the slope da/ din L approaches 
1/-7T (gray line) consistent with weak antilocalization. Here 
goo = g for dotted lines and goo = for dashed lines. For all 
other a/3, g a p — g. (Inset) Raw data a vs. Z//£. 



holds generally for WTFs (cf. |App. Dj ). 

The numerical data supporting the phase diagram in 
Fig. [1] are shown in Figs. [5] and [3] At in = the conduc- 
tivity always flows to the symplectic metal, regardless of 
the strength of the disorder (cf. Fig. [2|. By rescaling the 
length (L/£ — > L/C), we can collapse all the data on a 
single curve demonstrating one-parameter scaling along 
the m = line. At large conductivity, the beta function 
/3(<t) = d(\n a) / d(\n L) approaches l/ira as predicted for 
weak antilocalization 

By varying m, it is possible to drive the system to an 



insulator, as shown in Fig. 3(a) At small m, the system 



remains a symplectic metal. At some critical m, a metal- 
insulator transition occurs and it ceases to conduct. For 
a fixed nonzero m such that the clean system is insulat- 
ing, disorder drives the system into a metallic phase at 
some critical disorder strength g c that depends on m, as 
demonstrated in Fig. |3(b)| In both these cases, at large 
conductivity the slope da/dhiL approaches 1/7T, indica- 
tive of weak antilocalization. 

Conditions for localization. Since a WTI is always con- 
ducting in the absence of mass, it is pertinent to discuss 
under what circumstances one expects a nonzero mass. 
The potential V yz (r) couples valleys centered at different 
momenta and thus requires short-range scatters. Fur- 
thermore a nonzero mass can arise only when the surface 
potential is commensurate with an even number of unit 
cells, such as in the case of cleaving the surface at a crys- 
tal plane [l2|, or when the WTI is grown on a lattice- 
matching substrate. As such, a nonzero mass would be 
marked by an enlargement of the unit cell and would ap- 
pear in a crystal diffraction experiment as a peak of order 
G^/2, where G„ is a reciprocal lattice vector character- 
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FIG. 3. (a) Metal-insulator transition as m is varied. Con- 
ductivity is plotted vs. system size for fixed /x£ = 1 and 
9xo = gyz = 2. For large m the system flows to an insulating 
state, while for small m the system is conducting. Among 
the conducting curves, the slope da/ din L approaches 1/V at 
large a. The data show a metal-insulator transition at a c 
consistent with the known value of 1.42e 2 //i [TJ. (b) Metal- 
insulator transition as disorder strength g a p = g is varied. 
The plot is a vs. L/£ for fixed m£ = 0.05 and \i — 0. Increas- 
ing disorder increases the conductivity, inducing a transition 
from an insulating phase to a metallic one at some critical g. 
The dashed line indicates a slope of 1/tt. These figures are 
consistent with the phase diagram in Fig. [1] 



izing the weak topological invariants of the WTI |20j . On 
the other hand, a period-doubling perturbation could in- 
dicate a valley-mixing term other than m (the other pos- 
sible terms being {V x q}, (V yx ) or (V yy )). It may also be 
possible to measure in via spin- and angle-resolved pho- 
tocmission spectroscopy (spin- resolved ARPES), by com- 
paring the spin-up and spin-down intensities at wavevec- 
tor G„/2. This proposal is motivated by the form of 
the potential T v a z , which differentiates the up and down 
spins. Localization may also occur due to lattice effects 
or higher order terms in the Hamiltonian [l7, 21, 22 1. 

In the case where the WTI consists of an odd number 
of QSH layers, we argue that the mass must be iden- 
tically zero. Consider stacking n QSH layers, with each 
layer in the ai, a 2 plane, and the layers a 3 offset from one 
another. For simplicity, we impose a periodic boundary 
condition in the a 3 direction. The surface spectrum of 
a plane parallel to a 3 will have two Dirac cones, cen- 
tered on different time-reversal invariant momenta k a 
and kfo. such that (k& — k ) • a 3 = (G„/2) • a 3 = ir. 
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The second quantized kinetic Hamiltonian will be of the 
form \I>t (k - k a ) ■ cr* a + *£(k - k&) • cr* 6 . & and * 
are the creation and annihilation operators satisfying the 
boundary condition fy(r + n& 3 ) = ^(r). To cast this into 
the form of the effective Hamiltonian ((TJ), we perform 
the gauge transformation vE^ — !> ^E , A1 e 4kf ' ' r for each of the 
fermion species. The gauge transformation will, in gen- 
eral, change the boundary condition for the operators 
^ a and Notice that exp[i(k a - k 6 ) • (na 3 )] = (-1)", 
and hence for odd n the transformed operators will have 
differing boundary conditions: i.e., one periodic and one 
antipcriodic. The mass term coupling the fermion species 
together in the effective Hamiltonian must have antipcri- 
odic boundary conditions and, hence, averages to zero. 
Therefore, for an odd number of stacked QSH layers, m 
is zero and the surface (parallel to the stacking direction) 
always flows to a metallic phase. These results settle the 
question of which of the two possible flow diagrams con- 
sistent with the quasi- ID numerics in Ref. Il2| is actually 
realized. 

Two-parameter scaling. The existence of one- 
parameter scaling along the line m = suggests that 
there might be a two-parameter scaling collapse for the 
entire range of parameters when the mass is nonzero, 
analogous to the quantum Hall transition (in the A 
class) [Hl-IIil]. Figure 4(a) shows a possible flow for con- 
ductivity a and the (unknown) second scaling parameter 
j. The horizontal scale j distinguishes between the two 
topological phases, much in the same way as the Hall 
conductivity in the quantum Hall case. 

Even without a precise definition of j as an experi- 
mental quantity, we may still infer a number of proper- 
ties of the flow diagram, (i) For large conductivity a, 
/3(cr) is positive and a flows upward towards infinity, (ii) 
There are two insulating stable fixed points (crosses) at 
(a,j) — (0, ±oo ) and regions which flows toward them 
(shaded regions), (iii) Consequently, there must be un- 
stable fixed points (dots) at j = ±oo which mark a metal- 
insulator transition, (iv) Near j = 0, the system must 
flow to a metallic phase, as there should not be a direct 
phase transition between the two insulating phases. Fig- 
ure 4(a) gives the simplest flow diagram consistent with 



these requirements. 

The two-parameter scaling of (a, j) implies that er(L/£) 
cannot be collapsed onto a single scaling curve (as in 
Fig. [2]), but onto a family of curves parameterized by a 
single variable x. The scaling form is 



° = f(L/e 



(5) 



where all the microscopic parameters m, fi, g a p, £, etc., 
determine the conductivity only via the two functions x 
and f • 

In Fig. |4(b)| we present the accompanied data for our 
two-parameter scaling hypothesis, by collapsing a vs. 
L/£* onto a family of curves. For each curve, the pa- 
rameters \x and g a /3 were fixed while m is varied until 




FIG. 4. (a) Two-parameter flow which captures the QSH- 
metal-insulator transition in the All class. The scaling vari- 
ables are a and j, the latter of which separates the normal or 
QSH insulator phases, (b) Numerical data cr(L/£*), demon- 
strating that the conductivity curves may be collapsed onto 
each other. The gray line is the data at m = from Fig. [2] 
The raw data and parameters are given in |App. E) 



<y(L) overlaps with the existing set of curves. The data 
show reasonable agreement with the scaling form ([5]). 

Quantum transport at the surface of a weak topolog- 
ical insulator thus shows a scaling structure similar to 
that of the quantum Hall plateau transitions. It should 
be possible to interpret experiments on weak topological 
insulators in terms of the above Dirac model and possi- 
bly to control the parameter m by choosing a substrate 
whose lattice potential generates the massive perturba- 
tion. In addition, the electronic structure of thin films of 
STI's can be mapped to the two Dirac cone system stud- 
ied here, with the tunneling between the surfaces taking 
the role of the mass ( cf |App. C[ ) . Our results should mo- 
tivate the search for WTI candidate materials, of which 
there are few. It remains to be seen if the two-parameter 
flow is generic to all noninteracting disordered systems 
in the symplectic class. 

We gratefully acknowledge insightful discussions with 
Christopher Mudry, Shinsci Ryu, and Ashvin Vish- 
wanath. This work is supported by FENA (R.M.), DOE 
BES (J.H.B.), and NSF DMR-0804413 (J.E.M.). 
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A. Renormalization flow away from clean Dirac point 

Disorder average, either using the replica trick or a supersymmetric representation, gives arise to an interacting field 
theory with coupling constants given by the amplitudes g a p of the disorder correlator, as well as single particle 
potentials with amplitudes \i and m [26l . 127^ . The clean Dirac point g a p = m = /i = is a fixed point of this theory. 
Under renormalization the coupling constants flow away from the clean Dirac point according to the equations: 

dgoo 



goo (goo + 9x0 H~ Qyx 

+ gyy + 9yz + 9zo) + 9yz(9yx + 9yy), 

gxo(goo + gxo - g V x - g yv - g yz - g z o) + gzo(g y x + g yy ), 



dlnL 
. dg x0 
d\nL 

dg yx dg. 



dlnL ~9xog z o + g yz goo~: dlnL , 

, dgyz 

d\nL 
. dg z0 

dlnL 



g yz (-goo + gxo + g yx + g yy - g yz + g z o) + goo(g y x + g yy ), ( 6 ) 
g z o{goo - g x o — g yx — g yy - g yz + g z o) + g x o(g y x + g yy ), 



_ M + 7^(300 + 9x0 + g y x + 9yy + 9yz +9 Z 0>, 



dlnL 2?r 
dm m 



(—900 + 9x0 + 9 y x + gyy — 9yz + <7zo)- 



dlaL 2tt 

These equations have been obtained in one loop following the standard procedure used for Dirac fermions 28-3l|. 
Since m initially flows away faster than the disorder couplings, this suggest the shape of the phase diagram around 
the clean point is as in Fig. 2 in the main text. Other interesting features of these equations are that both the 
random scalar potential 500 and random mass g yz tend to decrease the mass, and thereby enhance the conductivity. 
In contrast, all the other disorder structures tend to increase the mass. 

Apart from these observation and that the clean fixed point is unstable, the renormalization group equations ((B]) 
yield no more information on the metal-insulator transition. Therefore the numerical simulations reported on in the 
main text are crucial in order to answer such questions. 



B. Chirality of the Dirac cones 



In the main text we used a model in which the two Dirac cones have the same chirality: 

a x k x + (T y k y 



H = Hv D 



a x k x + a y k % 



(7) 



Here H refers to the kinetic portion of the Hamiltonian. As the momentum k rotates by 2ir, the spin also rotates by 
2tt in the same direction, hence both Dirac cones have chirality of +1. 

We can apply the unitary transformation U which flips only one of the Dirac cone's chirality. 



U 



1 



H' = UHqU^ = hv D 



a x k x + a v k y 



<r x k x — a y k y 



(8) 



It is important to note that the transformation does not alter the form of time-reversal, i.e. UQU^ = G = i<J v 1C, and 
that it shuffles the disorder potentials: 



(KcO) Vy X , Vyy, Vy Z , V Z Q , VqQ ) ~ > (Vy X , ~ V X Q , ~ V yZ , Vyy , V Z Q , Vqq) ■ 



(9) 



The mass for the new system is defined to be m = (V yy ) as t v <j v anticommutes with H' . This shows that our results 
are independent of the chirality of the Dirac cones. 



C. Mapping to strong topological insulators thin films 

A thin film of strong topological insulator consists of two Dirac cones, for which the two cones will couple to each 
other for a sufficiently thin sample. The low-energy Hamiltonian for the thin film (with inversion symmetry) [32 , 
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is 



Hi 



STI 



hvD& x k 
A 



A 

-hyper x k 



+ V(v), 



(10) 



where cr x k = cr x p y — a y p x , and A is the tunneling amplitude between the two surfaces. The Hamiltonian is written 
in the basis (tf, t^, b j, bj.), where t and b represents the excitations for the top and bottom surfaces respectively. 

The model for WTI may be mapped to i?sTi v i a the transformation 



#STi = UHW 



where U — exp [ij 



T G 



(Again note that UQU^ = = ia v K.) The mass corresponding to this system is m = (V^o) = A/hvn, which is a 
measure of the tunneling amplitude between the surfaces. Consequently, the physics described in the main text may 
also be realized experimentally by a thin film of strong topological insulator, where the film's thickness can be used 
to tune m. (The mass m may oscillate between being positive and negative as the thickness is varied, as shown in 
Ref. SI) 



D. Constructing the mass for an arbitrary even number of Dirac cones 



For the case of two Dirac cones, we have defined the mass m = (V yz ) with the following properties: (1) In absence of 
other potentials and disorder, the system gap is simply 2|m|, and (2) the limits m — > oo and m — >■ — oo correspond to 
the system being in the trivial and QSH insulating phases, respectively. Here we explicitly construct m for a system 
with four Dirac cones and give the procedure for finding m in the general case. 

With four Dirac cones, the Hamiltonian is 

H& = v Q T°(T-k + V{r). (12) 



We use the Pauli matrices r x - y - z and v x,v,z acting in valley space to span all the possible intervalley couplings. 
(t° = v° = 1.) The disorder potential decomposes 

^(r)=^^ 7 (rKrV, with a, / 8, 7 6 {0, x, y, z). (13) 

Of the twenty-eight disorder structures w q t /3 ct 7 compatible with time-reversal Q = ia v lC : only six anticommutes with 
the kinetic Hamiltonian cr • k. These six are given by (a, /3, 7) = (0, y, z), (y, x, z), (y, z, z), (y, 0, z), (x, y, z), (z, y, z). 
The first three anticommute with each other as do the last three, while all of the first three commute with any of the 
last three. The mass is thus given by: 

™ (4) - (V 0yz ) 2 + (V yxz f + (V yzz ) 2 - yJ(V y0z ) 2 + (V xyz ) 2 + (V zyz ) 2 . (14) 

It can be shown that this definition satisfies (1) and (2) above. 

In the general case with 2n Dirac cones, there will be n(2n — 1) linear independent disorder structures which 
opens a gap in the system, of the form V ® <J Z , where V is an 2n x 2n matrix acting in valley space. V must be 
an antisymmetric pure-imaginary matrix, i.e. an element of the Lie algebra so(2n) (in the canonical representation). 
Seeing that V ® a z anticommutes with cr ■ k, the spectrum at k = is the set of eigenvalues of V. The mass is defined 
as follows: \m\ is the smallest non-negative eigenvalue of V while the sign of m is that of Pi(iV). 
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E. Raw data for two-parameter scaling 
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0.7 0.02 0.036 0.05 
0.5 1 0.194 0.25 0.299 
0.7 5 0.21 0.3 0.4 

1 5 0.36 0.485 0.618 
1 2 0.65 0.713 0.765 
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(Left) Numerical data er(Z//f *), demonstrating that the conductivity curves may be collapsed on to each other. The gray line is 
the data at m — from Fig. 2 of the main text. (Right) Raw data a(L/^). For each data set (denoted by color), the parameters 
£, n, g a p are fixed. The black data set has g x o — g y , = g, while for other colors, g a p = g for all a/3. Within each data set, m 
is varied until <r(L/£*) fit on top of one another. The numerical values are given as mi, mj, 7713. (Increasing m decreases the 
conductivity, hence the top curve of a data set has the smallest m.) 



F. Conductance Distribution 



0.8r (b) 




0.8r ( c ) 



9/{e 2 /h) 




(f) 



histogram 


m£ 


L/t 


samples 


(a) 


0.65 


160 


3456 


(b) 


0.67 


80 


9600 


(c) 


0.68 


160 


3456 


(d) 


0.7 


160 


1920 


(e) 


0.72 


80 


9600 



= 1, g*o = g V z = 2. 



(a)-(e) The conducta nce d istribution P as a function of the conductance g for a square geometry (W = L). The parameters 
used are those of Fig. 3(a) with varing m, given in the table (f). Comparing with the critical conductance distribution in Ref. [341 
computed for network models, we see that the metal-insulator transition occurs between (b) and (c) (0.67 < m£ < 0.68). This 
confirms that the 2D Anderson transition in the symplectic class is universal among both network models and continuous 
models with smooth disorder. 
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